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Abstract 

Measurements are presented of single and double-differential dijet cross sections in 
diffractive photoproduction based on a data sample with an integrated luminosity of 
47 pb _1 . The events are of the type ep — > eXY, where the hadronic system X con- 
tains at least two jets and is separated by a large rapidity gap from the system Y, which 
consists of a leading proton or low-mass proton excitation. The dijet cross sections are 
compared with QCD calculations at next-to-leading order and with a Monte Carlo model 
based on leading order matrix elements with parton showers. The measured cross sections 
are smaller than those obtained from the next-to-leading order calculations by a factor of 
about 0.6. This suppression factor has no significant dependence on the fraction x 7 of the 
photon four-momentum entering the hard subprocess. Ratios of the diffractive to the inclu- 
sive dijet cross sections are measured for the first time and are compared with Monte Carlo 
models. 
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1 Introduction 



Quantum Chromodynamics (QCD), as the gauge field theory of the strong interaction, reliably 
predicts scattering cross sections involving short distance partonic interactions. However, the 
vast majority of hadron-hadron scatterings take place through long-distance strong interactions, 
where no hard scales are present and perturbative QCD calculations are not possible. Promi- 
nent among these soft interactions are diffractive processes, in which the interacting hadrons 
remain intact or dissociate into low mass hadronic systems via an exchange which has vacuum 
quantum numbers, often referred to as a 'pomeron' Q]|. Following the observation of diffractive 
pp collisions in which a hard scale is provided by high transverse momentum jets [2], it has 
become possible to describe some classes of diffractive processes in terms of partonic interac- 
tions [3J. More recently, diffractive deep-inelastic scattering (DDIS) processes at HERA 01, of 
the type ep — > eXp, have been studied in detail and have led to a new level of understanding of 
the properties and structure of the diffractive exchange. Developing this microscopic descrip- 
tion of diffraction in terms of QCD and parton dynamics is a step towards a more complete 
understanding of the strong interaction. 

In the framework of a collinear factorisation theorem Q for hard scattering in semi-inclusive 
processes such as DDIS, diffractive parton distribution functions (DPDFs) may be defined. The 
DPDFs have similar properties to the standard parton distribution functions (PDFs) of the pro- 
ton, but with the constraint that there be a leading proton present in the final state. This condition 
may be satisfied equivalently by the experimental signatures of either a leading proton [6l[7]| or 
the presence of a large gap in the rapidity distribution of final state hadrons, separating an unob- 
served outgoing proton from the remainder of the hadronic final state In various extrac- 
tions using recent DDIS data lf8]4TTll. the DPDFs have been found to be dominated by gluons. 
To good approximation they exhibit a 'proton vertex factorisation' property, whereby they vary 
only in normalisation with the four-momentum of the final state proton, the normalisation being 
well modelled using Regge phenomenology lfT2l . 

Given a knowledge of the DPDFs, perturbative QCD calculations are applicable to other 
DDIS observables. Such calculations have been successful in the prediction of jet [|9l [T3l[T4l 
and heavy quark [15 ] production in DDIS at HERA. In both cases, next-to-leading order (NLO) 
QCD predictions using the DPDFs from [8] describe the measured cross sections well. How- 
ever, as has long been anticipated |[5lfT6ll. DPDF-based predictions for hard diffractive processes 
such as dijet production in pp scattering fail by around an order of magnitude to describe the 
data II 17111 81 . This factorisation breaking is generally attributed to absorptive corrections, corre- 
sponding to the destruction of the rapidity gap due to multiple interactions within a single event. 
Such effects are possible in pp scattering, where a beam remnant is present, in contrast to the 
electron scattering case in DDIS at HERA. A diversity of models of the absorptive corrections 
has developed HHI201, several of which reproduce the approximate 10% 'rapidity gap survival 
probability' observed in single diffractive pp scattering. 

The issues of DPDF applicability and rapidity gap survival can be studied in ep collisions at 
HERA in hard diffractive 'photoproduction', where the virtuality Q 2 of the exchange photon is 
close to zero. Under these circumstances, the photon can develop an effective partonic structure 
via 7 — > qq fluctuations and further subsequent splittings [21]. In a leading order picture, 
there are thus two classes of hard photoproduction: 'direct' interactions, where the photon 
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enters the hard scatter as a structureless object and 'resolved' interactions, where the photon 
interacts via its partonic structure and only a fraction x 7 of its four-momentum participates in the 
hard subprocess. Resolved photoproduction interactions can be further divided into a 'hadron- 
like' contribution and an 'anomalous' or 'point-like' contribution, the latter arising from the 
inhomogeneous term in the DGLAP equation for the photon ll22l . Interactions involving the 
hadron-like component resemble hadron-hadron scattering to a large extent and are therefore 
widely expected to exhibit gap destruction effects. The rapidity gap survival probability for 
these hadron-like processes has been estimated in a phenomenological model to be 0.34 H23 11241 . 
The point-like contribution to photon structure is expected to be subject to smaller absorptive 
corrections than the hadron-like part [|25l . In a recent model [|24l a survival probability of 
around 0.7 — 0.8 was obtained for diffractive dijet photoproduction, depending slightly on the 
jet transverse energies (Ej, ). 

Previous HI measurements of diffractive dijet photoproduction lfi~3ll26l have found cross 
sections to be smaller than NLO theoretical predictions, suggesting rapidity gap survival prob- 
abilities of around 0.5 with little dependence on x r A recent ZEUS measurement at somewhat 
larger E^ 11271 yielded a larger survival probability, compatible with unity. It has been proposed 
that this apparent discrepancy may be resolved if the rapidity gap survival probability depends 
on the scale of the hard interaction, an idea which is supported to some extent by data [|26] - l28l . 
Neither HI nor ZEUS data provide any evidence for the expected a; 7 dependence of the rapidity 
gap survival probability. 

A measurement of the ratio of the diffractive to the inclusive dijet photoproduction cross 
sections was proposed in H23[|25[|29l as a means of evaluating the gap survival probability. This 
ratio is expected to be relatively insensitive to the model of the photon parton densities and also 
offers cancellations of experimental systematics and higher order QCD corrections. A similar 
ratio was measured by the CDF collaboration ifPTl as a means of extracting effective pp DPDFs 
for comparison with HERA predictions and assessment of gap survival probabilities. 

This paper reports diffractive dijet photoproduction cross section measurements based on 
a positron-proton scattering data sample with luminosity about a factor three larger than that 
previously published by HI [|26l . The larger sample makes double-differential measurements 
possible, giving greater detail on the dynamics of gap survival and allowing studies of the 
correlations between the kinematic variables. The hypothesis of an E^ dependent rapidity gap 
survival probability is tested. The ratio of the diffractive to the inclusive dijet photoproduction 
cross sections is also extracted for the first time. 



2 Kinematic Variables 

Figures \T\ (a) and (b) show leading order examples of direct and resolved diffractive dijet pro- 
duction. Denoting the four-vectors of the incoming positron, the incoming proton and the ex- 
changed photon as k, P and q, respectively, the standard DIS kinematics can be described in 
terms of the invariants 



s = (k + P) 2 



Q 2 = ~q 2 



q ■ P 



Q 2 



(1) 



y = 



k-P 



x = 



2q- P 
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Here, s is the square of the total centre of mass energy of the collision, Q 2 is the photon virtu- 
ally, y is the scattered positron inelasticity and x is the fraction of the proton four-momentum 
carried by the quark coupling to the photon. These variables are related through Q 2 = s xy and 
to the invariant mass W of the photon-proton system by 

W = V(q + P) 2 « y/ys " Q 2 • (2) 




Figure 1 : Leading order diagrams for diffractive dijet photoproduction at HERA. Diagrams (a) 
and (b) are examples of direct and resolved photon interactions, respectively. 



Defining P x and P Y to be the four-vectors of the two distinct final state systems, where Y 
may be either a proton or a low mass proton excitation, the diffractive kinematics are described 
by the variables 

M\^P 2 X M 2 =P 2 t^(P-Py) 2 X]p EE — ~ . (3) 

Here, M x and M Y are the invariant masses of the systems X and Y, t is the squared four- 
momentum transfer at the proton vertex and xp is the fraction of the longitudinal momentum 
of the proton transferred to the system X. 

With u and v being the four-momenta of the particles entering the hard scatter from the 
photon and proton, respectively (u = q in the direct photon case), the invariant mass M i2 of the 
dijet system and the fractional photon (x 7 ) and pomeron (zp) longitudinal momenta entering 
the hard subprocess can be expressed as 

n ^7 P-u q-v 

M 12 EE y/{u + v) 2 x ^—q ZP = q-{P-Py) - () 

3 Theory and Models 

3.1 Diffractive Dijet Photoproduction in the Factorisation Approach 

Dijet electroproduction cross sections for Q 2 — > can be calculated in a fixed order QCD 
approach, assuming QCD collinear factorisation and neglecting any rapidity gap destruction 
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effects, as a convolution of partonic cross sections, photon PDFs and DPDFs, according to 

da(ep — > e + 2 jets + X' + Y) = / dy f 7 / e (y) / dx 1 fj/ 1 (x 1 , fip) ® 

J dt J dxp J dz F da(ij -> 2 jets) ff/ p {zp, fi F , x F ,t) . (5) 

Here, the hadronic system X' corresponds to the remainder of the system X after removing 
the two jets. The sum runs over all partons % and j that contribute, / 7 / e is the equivalent flux 
of photons that emerge from the incoming lepton 11301 and / 3 y 7 are the photon PDFs (/ 3 y 7 = 
5(1 — Xj) in the direct photon case). The hard partonic cross sections are denoted a, fP, are the 
DPDFs of the proton and fip is the factorisation scale. 

In this analysis, the GRV HO OTTl parton densities are used to describe the structure of the 
resolved photon. The HI 2006 Fit B set is used for the DPDFs, obtained by the HI collaboration 
in fits to inclusive DDIS data [8]. In the poorly constrained large Zp region, previous DDIS final 
state data flUIElllIl have shown a clear preference for these DPDFs over the Fit A set from 0. 
Further DPDF sets from HI (HI 2007 Fit Jets) [9] and ZEUS (ZEUS DPDF SJ) f7\\M in 
which dijet data from DDIS are used to improve the sensitivity to the large zp region, are also 
considered. All of these DPDF sets assume proton vertex factorisation, such that 

fi/ P (zp, Hf, xip, t) = f F / p (x P , t)fi/]p(zjp, fip) . (6) 

In the interpretation illustrated in figure [H fp/ p (x]p, t) may be considered as a pomeron flux, 
parameterised using Regge phenomenology in the DPDF sets used here. The fi/p factor then 
represents the parton densities of the pomeron. At relatively large values of xp > 0.01 a small 
contribution from a sub-leading meson exchange is required in the DPDF fits. This is taken 
into account by adding a second term of the same form as equation [6l but with different parton 
densities and a flux factor which is suppressed as xp — > 0. 



3.2 Next-to-leading Order Parton Level QCD Calculations 

The dijet electroproduction cross sections (equation [5]) are calculated at NLO of QCD using 
the program of Frixione et al. (FR) Ii32l adapted for diffractive applications as described in 
H26H33H . In the FR program, the renormalisation and factorisation scales are set to be equal and 
both are taken here from the leading jet transverse energy, i.e fin = fip — . The NLO 
calculations are performed with the number of flavours fixed to 5 and the QCD scale parameter 
set to A 5 = 0.228 GeV, corresponding to a 2-loop a s (M z ) = 0.118. The sensitivity of the 
calculated cross sections to the chosen fi R and fip values is studied by varying both scales 
simultaneously by factors of 0.5 and 2. The NLO calculations were cross-checked with the 
program written by Klasen and Kramer [25], which yields consistent results Il28ll3"4ll . 



3.3 Monte Carlo Simulations 

3.3.1 Corrections to the Data 

A Monte Carlo (MC) simulation is used to correct the data for detector effects in obtaining cross 
sections at the level of stable hadrons. All MC samples are passed through a detailed simulation 
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of the HI detector response based on the GEANT program 11351 and are subjected to the same 
reconstruction and analysis algorithms as are used for the data. 

Diffractive dijet photoproduction events are generated using the RAPGAP MC genera- 
tor 11361 in the range Q 2 < 0.01 GeV 2 with the minimum transverse momentum of the partons 
entering the hard subprocess set to p™ m = 2 GeV. RAPGAP is based on leading order QCD 
matrix elements with DGLAP parton showers. The HI 2006 fit B DPDFs and the GRV-G LO 
photon parton densities OT1 are used at a factorisation scale given by Ej, * . A sub-leading me- 
son exchange is included in the DPDF simulation, though its contribution is smaller than 5% in 
the kinematic range covered here. 

The selection of diffractive events using the rapidity gap method (section |4~2T) yields a sam- 
ple which is dominated by elastically scattered protons, but which also contains an admixture 
of events in which the proton dissociates to low M Y states. The measurement is corrected to the 
region M Y < 1.6 GeV and \t\ < 1 GeV 2 (section |43T ), using MC samples generated using the 
DIFFVM 11371 program, with and without proton dissociation, following the method described 
in©. 

In order to estimate the small (~ 2%) contribution to the data sample from non-diffractive 
events which pass the diffractive selection, the PYTHIA MC generator [381 is used in photo- 
production mode (Q 2 < 0.01 GeV 2 ). The PYTHIA model is also used to correct the inclusive 
photoproduction dijet measurement for detector effects when extracting the ratio of diffractive 
to inclusive cross sections. To estimate the model uncertainties on the corrections to the inclu- 
sive dijet cross section, a further MC sample is obtained using the HERWIG generator [|39ll . 
Multiple interaction models are included in both PYTHIA and HERWIG as described in sec- 
tion 13.3.21 With these settings, PYTHIA provides a good description of the shapes of the un- 
corrected inclusive dijet distributions, with a normalisation slightly larger than that of the data. 
The HERWIG MC underestimates the normalisation of the cross section by about a factor of 
two, but gives an acceptable description of the shapes of the measured distributions. 



3.3.2 Corrections to Theoretical Models 

For comparison with the diffractive measurements, it is necessary to convert the calculated NLO 
parton-level cross sections to the level of stable hadrons by evaluating effects due to initial and 
final state parton showering, fragmentation, hadronisation and the influence of beam remnants. 
The RAPGAP MC is used to compute the required 'hadronisation correction' factors to the 
diffractive dijet calculations. These factors are defined for each measured data point by 

^.hadron 

1 + Shadr. = parton ■ ( 7 ) 
^dijet 

They reduce the predicted cross sections by typically 15% and are given for each data point in 
tables [2] and HI The shape of the x y distribution is most strongly affected by the hadronisation 
corrections, the main effect being the migration of some direct photon interactions, for which the 
cross section is large, to lower x 7 values, substantially increasing the prediction in the interval 
0.6 < s 7 < 0.8. 
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Resolved photon interactions in inclusive dijet production [401 are poorly described by NLO 
calculations unless a model of multiple interactions (MI) within a single event is included in ad- 
dition to the low transverse momentum beam remnant, QCD radiation and hadronisation con- 
tributions which are simulated in standard MC models. The PYTHIA MC generator is used to 
investigate the influence of Mis on inclusive dijet cross sections and hence on the diffractive-to- 
inclusive ratios. Several tunes of the PYTHIA model for multiple hard parton-parton scatterings 
are available. In the version used here (tune A in [38]), both the proton and the resolved photon 
have double-Gaussian matter distributions. The minimum transverse momentum down to which 
secondary partonic scatterings are calculated depends on the impact parameter of the collision 
and is governed by a regularisation scale p± = 1.2 GeV |[38l . 

Multiple hard partonic interactions are not simulated in the HERWIG MC. Instead, an em- 
pirically motivated soft multiple interaction model [[391 is used. The probability of such activity 
in a given event is set to 0.25. 



4 Experimental Procedure 
4.1 The HI Detector 

A detailed description of the HI detector can be found elsewhere iflTTl . Here, a brief account 
is given of the detector components most relevant to the present analysis. The HI coordinate 
system is defined such that the origin is at the nominal ep interaction point and the polar angle 
9 = 0° and the positive z axis correspond to the direction of the outgoing proton beam. The 
region 9 < 90°, which has positive pseudorapidity 77 = — lntan6'/2, is referred to as the 
'forward' hemisphere. 

The ep interaction point in HI is surrounded by a central tracking region, which includes 
silicon strip detectors as well as two large concentric drift chambers. These chambers cover 
a pseudorapidity region of —1.5 < 77 < 1.5 and have a transverse momentum resolution of 
<j(Pt)/Pt — 0.006PT/GeV © 0.02. They also provide triggering information. The central 
tracking detectors are surrounded by a finely segmented liquid argon (LAr) sampling calorime- 
ter covering —1.5 < 77 < 3.4. Its resolution is a / E = 0.11/ yE~/ GeV © 0.01 for electrons and 
photons and a / E = 0.50/ WE] GeV © 0.02 for hadrons, as measured in test beams R2ll . The 
central tracker and LAr calorimeter are placed inside a large superconducting solenoid, which 
produces a uniform magnetic field of 1.16 T. The backward region — 4 < 77 < — 1.4 is covered 
by a lead-scintillating fibre calorimeter (SpaCal) with electromagnetic and hadronic sections. 

Information from the central tracker and the LAr and SpaCal calorimeters is combined using 
an energy flow algorithm to obtain the hadronic final state (HFS) [|43l . The hadronic energy 
scale is known to 3% for this analysis |[34l . 

Photoproduction events are selected by tagging positrons scattered through very small an- 
gles, corresponding to quasi-real photon emission, using a crystal Cerenkov calorimeter at 
z = —33 m (electron tagger). The luminosity is measured via the Bethe-Heitler bremsstrahlung 
process ep — » e/ry, the final state photon being detected in another crystal calorimeter at 
z = -103 m. 
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A set of drift chambers around z = 6.5 m comprise the forward muon detector (FMD). 
The proton remnant tagger (PRT) is a set of scintillators surrounding the beam pipe at z = 
26 m. These detectors, used together with the most forward part of the LAr, are efficient in the 
identification of very forward energy flow and are used to select events with large rapidity gaps 
near to the outgoing proton direction. 



4.2 Event Selection, Kinematic and Jet Reconstruction 

The analysis is based on a sample of integrated luminosity 47 pb _1 , collected by HI in 1999 
and 2000 with proton and positron beam energies of 920 GeV and 27.5 GeV, respectively. The 
events are triggered on the basis of a scattered positron signal in the electron tagger and at least 
three high transverse momentum tracks in the drift chambers of the central tracker. 

The event inelasticity y and hence the invariant mass W of the photon-proton system are 
reconstructed using the scattered positron energy E e > measured in the electron tagger according 
to 

y = l -^L w = ^y, ( 8) 

e 

where E b e eam is the positron beam energy. The geometric acceptance of the electron tagger 
limits the measurement to Q 2 < 0.01 GeV 2 and intermediate values of y. 

The reconstructed hadronic final state objects (section 14.11) are subjected to the k T longi- 
tudinally invariant jet algorithm [44J, applied in the laboratory frame with parameters R = 1 
and E 3 ^ t,mm = 2.5 GeV. To facilitate comparisons with the NLO calculations, different cuts 
are placed on the transverse energies Erj? and E^ t2 of the leading and next-to-leading jets, 
respectively. As well as these variables, the jet properties are studied in terms of the variables 

| Ar/ jcts | = |?? jotl - rf ct2 \ <r/ jots ) = - (rf ctl + rf ct2 ) , (9) 

obtained from the laboratory frame pseudorapidities of the jet axes. With J^ 1 ' and denoting 
the four-momenta of the two jets, hadron level estimators of the dijet invariant mass and of x 1 
are obtained from 

M 12 = V2 JW ■ J® ay = ^ jets{ f' ~ , (10) 

where the sums labelled 'HFS' and 'jets' run over all hadronic final state objects and those 
included in the jets, respectively. 

The diffractive event selection is based on the presence of a large forward rapidity gap. The 
pseudorapidity of the most forward cluster in the LAr calorimeter with energy above 400 MeV is 
required to satisfy r] max < 3.2. The activity in the PRT and the FMD is required not to exceed 
that typical of noise levels as obtained from randomly triggered events. These requirements 
ensure that the analysed sample is dominated by elastically scattered protons at small \t\, with 
a small admixture of events with leading neutrons and low My baryon excitations, collectively 
referred to here as 'proton dissociation' contributions. 



10 



The diffractive kinematics are reconstructed using 



^HFs{Ei + P Z) i) 

x *> = M x = ^fx^ } (11) 

where E h ^ am is the proton beam energy. A cut on x$> is applied to ensure good containment of 
the system X and to suppress sub-leading exchange contributions. A hadron level estimator for 
the momentum fraction z$> is obtained using 

jets ^jetsy-Ej -\- P z ,i) f\0\ 

** ~ ZHFs(Ei + P z>i ) ' ( ] 



The kinematic range in which the diffractive dijet measurement is performed is specified in 
tabled] The inclusive measurement phase space is defined by these conditions, with the require- 
ments relaxed on the diffractive variables xjp, z S p S , M Y and t. Except where measurements are 
made explicitly as a function of z 3 p S , the z 3 p S > 0.8 region is excluded from the diffractive anal- 
ysis. This improves the reliability of the comparison between data and theoretical predictions, 
since the DPDF sets used are not valid at the largest Zp values. After applying all selection 
criteria, about 3600 out of roughly 200000 inclusive dijet photoproduction events are used in 
the diffractive analysis. A more detailed description of the analysis can be found in ll34~l . 



Diffractive and Inclusive Measurements 


Q 2 < 0.01 GeV 2 
E^ 1 > 5 GeV 
-1 < if ctl < 2 


0.3 < y < 0.65 
E^ 2 > 4 GeV 
-1 < V ct2 < 2 


Diffractive Measurement 


x F < 0.03 
My < 1.6 GeV 


4 ts < 0.8 
\t\ < 1 GeV 2 



Table 1: Kinematic ranges of the diffractive and inclusive measurements. 



4.3 Cross Section Measurement 



The diffractive differential cross section is measured in each bin i of a variable x using the 
formula 

/da\ Nf a ta j^ 19 - Nf C ' bgd 1 

\dx) i Ai A? £ ' CP diss ' 

Here, Nf ata is the raw number of reconstructed events passing the selection criteria listed in 
section 14.21 and e l [ 19 is the trigger efficiency, obtained by reference to an independently trig- 
gered sample and parameterised as a function of the multiplicity of charged particle tracks. The 
trigger efficiency averaged over the full measurement range is 0.86. The non-diffractive back- 
ground contribution obtained from the PYTHIA MC simulation is denoted N* 1 ' b9d and does 
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not extend beyond the few percent level for any of the measured data points. The factor A,i 
corrects the measurement for detector effects, including migrations between bins, to the level 
of stable hadrons. It is calculated from the RAPGAP MC and has an average value of 0.31 for 
the diffractive analysis, most of the losses being due to the limited electron tagger acceptance 
of 0.40 integrated over the measured y range. The bin width is denoted Af , C is the luminosity 
of the data sample and C pdiss = 0.94 ± 0.07, evaluated using the DIFFVM MC, corrects the 
measurement to the chosen range of M Y and t (table [[])• In the inclusive analysis the cross 
section is obtained analogously to equation [13] except for the N MC ' bgd and C pdlss terms, which 
are not relevant. 

4.4 Systematic Uncertainties 

Uncertainties are evaluated for all significant sources of possible systematic bias. These sources 
are summarised for the diffractive analysis below, together with their corresponding influences 
on the total diffractive cross section. 

Energy Scale: The energy scale of the HFS measurement is tested using the momentum bal- 
ance constraint between the precisely reconstructed positron and the HFS in neutral current DIS 
events. Dedicated data and MC samples are analysed and found to agree to better than 3%. The 
effect of a relative 3% change in the energy of the HFS between the data and the MC is a 9.6% 
shift to the total diffractive dijet cross section. This arises mainly from changes in the migration 
corrections across the minimum Ej? values and the maximum xp value of the measurement. 
The energy scale uncertainties are thus highly correlated between the bins of the differential 
cross section measurements. 

Large Rapidity Gap Selection: A fraction of the events in the kinematic range of the analysis 
(table Q]) give rise to hadronic activity in the forward detectors or at pseudorapidities beyond 
those allowed by the i] max cut in the LAr calorimeter. Corrections for this inefficiency of the 
large rapidity gap selection are made using the RAPGAP MC simulation. The uncertainties in 
the correction factors are assessed through a study of forward energy flow in a sample of dijet 
photoproduction events with leading protons tagged in the HI Forward Proton Spectrometer H. 
RAPGAP is found to describe these migrations to within 10% 112611451 . which translates into a 
2.9% uncertainty on the measured total cross section and uncertainties which are correlated 
between bins of the differential distributions. 

Proton Dissociation: The model dependence uncertainty on the proton dissociation correc- 
tion factor (CP dlss in equation [T3~l) is obtained by varying the elastic and proton dissociation 
cross sections and the proton dissociation My and t dependences in the DIFFVM MC samples, 
following [8J. The largest effect arises from varying the ratio of the proton-elastic to the proton- 
dissociative cross sections between 0.5 and 2. The resulting uncertainty on the measured cross 
section is 7%. 

Model Dependence: The influence of the model assumptions on the acceptance and bin mi- 
gration corrections (Ai in equation [T3l). is determined in the diffractive analysis by varying the 
kinematic distributions in the RAPGAP simulation within the limits allowed by maintaining 
an acceptable description of the uncorrected data. The following variations are implemented 
by reweighting each MC event according to the value of generator level kinematic variables, 
leading to the quoted systematic uncertainties on the total cross section. 
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• The x P distribution is reweighted by x^>° 2 , leading to a 6.0% uncertainty. 

• The z*p S distribution is reweighted by z^ s ±0 ' 3 , leading to a 4.8% uncertainty. 

• The a; J 7 ets distribution is reweighted by x J 7 ets ±0,3 , leading to a 0.6% uncertainty. 

• The distribution is reweighted by E^ 1 ±0A , leading to a 0.8% uncertainty. 

• The t distribution is reweighted by e ±2t , leading to a 4.4% uncertainty. 

• The y distribution is reweighted by y ±0 3 , leading to a 0.2% uncertainty. 

Electron Tagger Acceptance: A dedicated procedure external to this analysis is used to ob- 
tain the electron tagger acceptance ll4~6l . The integrated acceptance over the full y range is 
known to 5%, which affects the cross section normalisation. 

Trigger Efficiency: The procedure for parameterising the trigger efficiency (ef 19 in equa- 
tion [13]) leads to a 5% uncertainty. This covers the observed deviations of the parameterisation 
from the measured efficiencies as a function of all variables relevant to the analysis. This un- 
certainty is treated as being uncorrelated between data points. 

Luminosity: The measurement of the integrated luminosity has an uncertainty of 1.5%. This 
translates directly into a 1.5% normalisation uncertainty on the measured cross sections. 

Non-diffractive background: A 50% normalisation variation is applied to the non-diffractive 
background contribution given by the PYTHIA MC model (N MC > b 9 d j n equation[[3]). The effect 
of this change is correlated between the data points and leads to a 1% uncertainty on the total 
cross section. 

Forward Detector Noise: Fluctuations in the FMD noise, leading to losses in the large ra- 
pidity gap event selection, are evaluated for each run using randomly triggered events. The 
standard deviation in the run-by-run distribution of the correction factors is used to derive a 
0.5% normalisation uncertainty on the measured cross sections. Noise in the PRT detector is 
negligible. 

A similar procedure is followed to evaluate the systematic uncertainties in the inclusive di- 
jet analysis. The uncertainties associated with the large rapidity gap selection and the model 
dependence are no longer relevant. Instead, comparisons between the PYTHIA and HERWIG 
MCs are used determine a 2% model dependence uncertainty on the acceptance correction when 
integrated over the full phase space studied. The inclusive cross section systematics are domi- 
nated by a contribution at the 10% level from the HFS energy scale uncertainty. However, when 
forming the ratio of diffractive to inclusive cross sections, this error source cancels to good ap- 
proximation, the residual uncertainty being less than 1%. The largest remaining contribution to 
the systematic uncertainty on the cross section ratio arises from the model dependence. 

The total systematic uncertainty on each data point is formed by adding the individual con- 
tributions in quadrature. In the figures and tables that follow, the systematic uncertainties are 
separated into two categories: those which are uncorrelated between data points (the model de- 
pendence and trigger efficiency) and those which lead to correlations between data points (all 
other sources). 
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5 Results 



5.1 Diffractive Dijet Cross Sections 

Cross sections are measured integrated over the full kinematic range specified in tableCQand also 
single- and double-differentially as a function of a variety of variables which are sensitive to the 
overall event structure, the hard subprocess and the presence of remnants of the virtual photon 
and the diffractive exchange. The measured differential cross sections, which correspond to av- 
erages over the specified measurement intervals, are given numerically in tables [2] and |4l where 
the experimental uncertainties and hadronisation corrections applied to the NLO calculations 
are also listed. Tables [3] and \5\ contain the ratios of the measurements to the NLO calculations, 
obtained using the FR framework ( section I3T21) and the HI 2006 Fit B DPDFs (referred to in the 
following as FR Fit B). 

5.1.1 Integrated Cross Section 

The total diffractive dijet positron-proton cross section integrated over the full measured kine- 
matic range is 

a^ ta = 295 ± 6 (stat.) ± 58 (syst.) pb . (14) 
The ratio of this result to the corresponding FR Fit B NLO prediction is 

PdatJvNLO = 0-58 ± 0.01 (stat.) ± 0.12 (syst.) ± 0.14 (scale) ± 0.09 (DPDF) , (15) 

where the statistical and systematic uncertainties originate from the measurement. The scale 
uncertainty corresponds to the effect of simultaneously varying the renormalisation and factori- 
sation scales from their central values, fi R ■ = \i F = E^ t , by a factor of two in either direction. 
This large (25%) scale uncertainty arises due to the relatively low Ej? range of this analysis, 
and is the limiting factor in the comparison between data and theory. The DPDF uncertainty 
is obtained using the method of 071, by propagating the eigenvector decomposition of the fit 
uncertainties. If the HI 2006 Fit B DPDFs are replaced by the HI 2007 Fit Jets DPDFs, the 
result is o^atJ ot NLO = 0.64, which is inside the quoted DPDF uncertainty. Using ZEUS DPDF 
SJ, a compatible result of o^atal at NLO = 0-^0 * s obtained. 

Adding all uncertainties in quadrature, the ratio result in equation Q3] implies at the 2a level 
that the NLO QCD calculation, neglecting any gap destruction effects, yields a larger diffractive 
dijet photoproduction cross section than that measured. It confirms the result of a previous HI 
analysis in a very similar kinematic range 11261 and is broadly as expected from theoretical 
calculations of rapidity gap survival probabilities H23H291 . 

5.1.2 Single-Differential Cross Sections 

Figure [2] shows the diffractive dijet cross section measured single-differentially in x J 7 ets , E^ * , 
logxjp, Zp S , (r] jets ), |Ar] jets |, W, M 12 and M x , in the phase space defined in table Q] In 
figureOd), the z 3 p S < 0.8 requirement is relaxed and the cross section measured for the region 
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z 3 jp S > 0.8 is shown without theoretical comparisons, since the DPDFs are not defined (see 
section |4T2l) . To allow a more detailed shape comparison between the data and the predictions, 
ratios of the measured differential cross sections to the FR Fit B calculations are plotted in 
figure [3l These ratios may be taken as measurements of the dependence of the rapidity gap 
survival probability on the kinematic variables. 

Figures [2] and [3] show that the suppression by around a factor of 0.6 of the data with respect 
to the FR Fit B NLO calculations has at most a weak dependence on the kinematic variables. 
Notably, within the uncertainties there is no dependence on x J 7 ets (figure [3£a)), in contrast to 
theoretical predictions for the rapidity gap survival probability H23 11291 . The largest dependence 
of the central values of the measured ratios on any of the variables appears in the cross section 
differential in (figure Ob)). Although not well established by the current data, this depen- 
dence is compatible with previous data [|26] - |28l . The dependence is investigated further 
in section [5TT31 

The measured cross sections in figure [2] are also compared with a prediction obtained using 
the RAPGAP MC generator (section |3 .3 .11) . which does not contain any model of rapidity gap 
destruction. The shapes of the measured cross sections are well described and the normalisation 
is only slightly lower than that of the data. However the scale uncertainty in this model is rather 
large and the same model undershoots diffractive dijet measurements in DIS (9l|26]|, where 
factorisation is expected to hold. 

In [|27l . the ZEUS collaboration presented an analysis of diffractive dijet photoproduction 
data with > 7.5 GeV, which is most readily compared with the second and third E^ 1 
intervals in figures (2b) and Ob). However, even for E^ % > 7.5 GeV, a direct comparison 
between HI and ZEUS data is not possible, since the ZEUS analysis covers a wider y range and 
cuts on the second jet at an E^ t2 value of 6.5 GeV, larger than the value used here. An indirect 
comparison can be made on the basis of ratios of the data to NLO theoretical calculations using 
the HI Fit B DPDFs. ZEUS obtains a result of around 0.9 for this ratio, which is compatible 
with the result for E^ 1 > 7.5 GeV obtained here, within the large combined uncertainties 

As discussed in detail in (HIED, the error bands on the DPDFs extracted from inclusive 
diffraction alone do not include uncertainties due to parton parameterisation choices and thus do 
not reflect the full uncertainties, particularly in the large z 3 p S region. To give a complementary 
indication of the possible range of variation, comparisons between the ratios obtained with the 
HI 2006 Fit B DPDFs, the HI 2007 Fit Jets DPDFs and ZEUS DPDF SJ fit are shown for a 
subset of variables (x J 7 ets , E^ n and z 3 p S ) in figurelU The ZEUS DPDFs lead to ratios which are 
uniformly 10 — 15% larger than those obtained with HI 2006 Fit B, with no strong dependence 
on any of the kinematic variables. The deviation of the HI 2007 Fit Jets result from the HI 
2006 Fit B result extends beyond the DPDF error band for z 3 p s > 0.6, which is correlated with 
a somewhat stronger dependence of the ratio of data to theory on E^ n and a slightly different 
shape at low x^ ts . 

According to the RAPGAP model, approximately half of the cross section in the kinematic 
range studied arises from each of the direct and resolved photon-induced contributions. The de- 
composition of photoproduction processes into direct and resolved interactions is not uniquely 

'in ||7]|48] a 13% difference between HI and ZEUS inclusive DDIS data is identified. This is within the 
combined normalisation uncertainties of the two experiments, which are largely due to proton dissociation. If the 
dijet photoproduction cross sections in the two experiments are normalised to the inclusive DDIS data ||7][8], the 
remaining differences in the common E^ 1 range are well within the experimental uncertainties alone. 
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defined beyond LO. When modelling rapidity gap survival probabilities in the following, the 
resolved photon contribution is defined to correspond exactly to that which is calculated using 
the photon structure function]^ Following the calculation using an absorptive model of a gap 
survival probability of 0.34 for the hadron-like component of resolved photoproduction [|23l . 
previous HI data [26j were compared in E51 with predictions in which the full resolved photon 
contribution was suppressed by this factor, the direct photon contribution being left unsup- 
pressed. In a later analysis ||29ll , this procedure was extended to NLO. The conclusions of these 
previous studies are confirmed in figures [5] and [6] through a similar comparison of the current 
data with NLO calculations in which the resolved photon contribution is globally suppressed 
by a factor of 0.34. The overall normalisation of this calculation is in good agreement with the 
data. However, the shapes of some of the differential distributions are not well reproduced. In 
particular, there is a variation by more than a factor of two in the ratio of data to theory as a 
function of x J 7 ets (figured). 

The distinction between point-like and hadron-like resolved photon interactions recently de- 
veloped in j[24| leads to a significantly weaker predicted suppression in the kinematic range of 
the current analysis. The data are compared with this refined 'KKMR' model under the approx- 
imation of completely neglecting hadron-like resolved photon contributions, which, according 
to the authors, become dominant only for x^ ts < 0.1 [50J, beyond the range of the current anal- 
ysis. The rapidity gap survival probabilities obtained in E4l for point-like photon interactions 
using the GRV HO photon PDFs are applied to all resolved photon interactions. Interactions 
involving quarks and gluons from the photon are thus suppressed by factors of 0.71 and 0.53, 
respectively. The quark-initiated contribution is dominant throughout the measured range, such 
that the rapidity gap survival probability in the model is approximately 0.7 for resolved photon 
interactions and 1 for direct photon interactions. 

Figure |5] shows the comparison between the measured single differential cross sections 
and the NLO QCD predictions, with the resolved photon contribution scaled according to the 
KKMR model. The corresponding ratios of data to theoretical predictions are shown in figure[6l 
The overall normalisation of the KKMR-based calculation is larger than that of the data, but is 
compatible within the large uncertainties. Many of the distributions studied are well described 
in shape (£^ etl , z } p S , |A?7 jets |, W and M l2 ). The data thus agree with the prediction [|29l that 
the Erj? dependence of the data/theory ratio flattens if the resolved photon contribution alone 
is suppressed. However, there remains a variation in the ratio of data to the KKMR model with 
x J 7 ets and to a lesser extent with (77 jets ), xjp and Mx- A comparison of figures [3] and [6] shows 
that the shapes of the differential cross sections are generally better described with a global 
suppression factor than with a survival probability applied to resolved photon interactions only. 

5.1.3 Double-Differential Cross Sections 

To study further the dynamics of rapidity gap suppression and their dependence on the nature of 
the photon interaction, cross sections are measured double differentially in two regions of x J 7 ets , 

2 In ||29|. an alternative procedure is introduced, whereby the part of the direct contribution which depends on 
the factorisation scale at the photon vertex is also suppressed, stabilising the dependence of the combined direct and 
resolved cross sections on this scale (49) . The difference between the rapidity gap survival probabilities obtained 
using the two methods (6% in fits |29l to previous HI data |26|) is small in comparison to other uncertainties. 
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which are enriched with either resolved (x J 7 cts < 0.75) or direct (x J 7 ets > 0.75) photon processes. 
Using the RAPGAP MC model with the GRV-G LO photon PDFs, the x J 7 ets < 0.75 region is 
estimated to contain 77% resolved photon interactions integrated over the measurement region 
(table Q]), with a 75% direct photon contribution for x^ ts > 0.75. 

In figures |7£a)-(c), measurements are presented of the double-differential dijet cross section 
d 2 cr/d_E'^ tl dx J 7 ets for three E^ ranges in the resolved and direct photon-enriched x J 7 ets intervals. 
The data are compared with the FR Fit B calculations and with the RAPGAP MC predictions. 
Due to kinematic constraints, the resolved-enriched cross section at low x J 7 ets falls most rapidly 

as increases. There is a suggestion that this E^ 1 dependence in the resolved-enriched 
region is stronger for the NLO QCD theory than for the data. In the direct-enriched high x^ ts 
region, the cross section falls more slowly with E^ t and the dependence in the data is similar to 
that predicted by the NLO calculation. These features are illustrated further in figures|7j(d)-(e), 
where the ratios of the data to the NLO theory from figures [7] (a)-(c) are presented as a function 
of E^ n in the resolved and direct photon-enriched x^ ts regions, respectively. 

The significance of the E^ dependence in the resolved-enriched region (figure [7J (d)) is 
evaluated through a x 2 test. All uncertainties are taken into account in this procedure, though 
the main contribution comes from the statistical and uncorrected systematic uncertainties on 
the data, the remaining uncertainties changing only the normalisation of the ratio to first ap- 
proximation. A test of the hypothesis that there is no E^ 1 dependence yields a x 2 value of 
1.36, with two degrees of freedom, corresponding to an E^ n variation at the 73% confidence 
level. The suppression of the data relative to the NLO prediction in the direct-enriched large 
^jets re gj on i S; within errors, independent of Ej? (figure [TJ (e)). Figures [TJ (d)-(e) thus indi- 
cate that any E^ n dependence of the data-to-theory ratio in figure Ob) is driven primarily by 
resolved photon interactions. An E^ n dependence of the gap survival probability is predicted 
in the KKMR model, due to variations in the size of the qq dipole produced by the point-like 
photon splitting, and hence in the absorptive correction. However the predicted effect is small 
(4% as E^ 1 changes from 5 GeV to 7.5 GeV). Figures dd)-(e) also indicate that when 
becomes large, the suppression in the direct region may be stronger than that in the resolved re- 
gion, which is not expected in any model. The large uncertainties permit statistical fluctuations 
in the data or small inadequacies in the theory as possible explanations. 

In figure [H the cross section is shown double differentially in z 3 p S and x^ ts . The measured 
cross section is compared with the NLO theory as a function of z^ s in two bins of x J 7 ets in fig- 
ures [8] (a)-(b) and the ratios of data to theoretical predictions are shown in figures [8] (c)-(d). The 
NLO calculations describe the measured shapes rather well, with no evidence for any variation 
of the suppression factor between any of the measurement ranges. The gap survival probability 
in a region where there are small or no remnants of either the photon or the diffractive exchange 
(highest z 3 p S bin in figure [8]i) is thus similar to that where both remnants are significant (lowest 
z 3 p S bin in figure[8]:). This remains the case when the HI 2007 Fit Jets DPDFs are used in place 
of HI 2006 Fit B. In both figures |7J and [8} the RAPGAP MC prediction gives a satisfactory de- 
scription of the shapes of the double differential cross sections, the normalisation being slightly 
lower than that of the data. 
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5.2 Ratios of Diffractive to Inclusive Cross Sections 



Measurements of ratios of diffractive to inclusive dijet photoproduction cross sections have been 
proposed Il2~3ll2~5ll2~9l as a further test of gap survival issues. Their potential advantages over 
straight-forward diffractive measurements lie in the partial cancellations of some experimental 
systematics and of theoretical uncertainties due to the photon structure and factorisation and 
renormalisation scale choices. The sensitivity to absorptive effects of diffractive-to-inclusive 
ratios is thus potentially superior to that of pure diffractive cross sections. For the ratio ex- 
traction presented here, inclusive dijet cross sections are measured using data collected in the 
same period as the diffractive sample. The experimental method and systematic error treatment 
for the inclusive case is described in section HI It is identical to the diffractive measurement 
method, with the exception of the large rapidity gap requirements. 

At the relatively low transverse energies studied in the present analysis, underlying event 
effects have a large influence on jet cross sections in inclusive photoproduction ll40ll . Here, the 
PYTHIA and HERWIG MC models are used to correct the inclusive data for detector effects, 
with MI included as described in section 13.3.21 The two models agree rather well on the cor- 
rections to be applied to the data. The average of the results with the two models is therefore 
used to calculate the corrections and the uncertainty. The latter is taken from the difference 
between the results with the two models and is relatively small (2% when integrated over the 
full measured range). 

The ratios of diffractive to inclusive single-differential dijet cross sections are given numer- 
ically in tableland are shown in figure [9] as a function of x J 7 ets , E^ 1 , (^ cts ), \ A^ ets \, M 12 and 
W. Due to the partial or complete cancellations of some error sources when forming the ratio, 
the correlated uncertainties are reduced compared with those for the diffractive distributions. 

Since they give adequate descriptions of the diffractive and inclusive data, respectively, the 
RAPGAP and PYTHIA MC models are used to assess the relative sensitivity of the diffractive- 
to-inclusive ratio to the gap survival and MI effects. With no MI effects included in the PYTHIA 
model, the description of the inclusive data is poor and the ratio of RAPGAP to PYTHIA 
exceeds the data by a factor of around 1.5. As expected, this factor becomes smaller as 
increases. However, the shape of the prediction also differs from that of the ratio data for most 
of the other variables studied, in particular x J 7 ets . 

The inclusion of the PYTHIA MI model changes the predicted inclusive cross sections, 
and hence the ratios, substantially. The ratio of RAPGAP to PYTHIA then gives an improved 
description of the shapes of the distributions. The MI effects alter the predicted ratio by a factor 
of 0.5 at low x J 7 ets , where the resolved photon remnant is most important. As expected, there 
is little effect in the direct photon-dominated large x J 7 ets region. The normalisation of the ratio 
of the models when MI are included is smaller than that of the data. This partially reflects the 
RAPGAP description of the diffractive data (figure [2]) and is partially due to an overshoot in the 
PYTHIA description of the inclusive data. 

The fractional reduction in the predicted inclusive cross section when MI are introduced in 
the PYTHIA model is comparable to the magnitude of the gap survival suppression factor in 
the diffractive data (section I5TTT) . The uncertainties in modelling the MI are large and difficult to 
quantify. The precision with which gap survival issues can be unfolded from MI complications 
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in the ratio of diffractive to inclusive data is correspondingly poor. Therefore no strong conclu- 
sions can be drawn with our current understanding of MI, despite the relatively good precision 
of the data. 

6 Summary 

Single and double-differential cross sections are measured for diffractive dijet photoproduction 
and are compared with predictions based on NLO QCD calculations using different sets of 
DPDFs. Ratios of the measured to the predicted differential cross sections are also studied. 

The total diffractive dijet cross section is overestimated by the NLO QCD theory by about 
a factor of two, which is consistent with previous HI measurements Il2~6l . The shapes of the 
single-differential cross sections are well described when the HI 2006 DPDF Fit B partons 
are used. A good overall description of the differential cross sections is obtained by applying 
a global suppression factor of 0.58 ± 0.21 to the NLO calculations. As in similar previous 
analyses Il2~6l - l2~8ll . there is a suggestion of a dependence of the rapidity gap survival probability 
on , though the significance of this effect is not large. 

If only the resolved photon contribution in the calculation is suppressed by a factor of 0.34, 
as predicted for hadron-like resolved photon interactions [23J, the overall normalisation of the 
NLO QCD prediction agrees well with the data. However, the description of the x J 7 cts distribu- 
tion, which best distinguishes direct from resolved photon interactions, becomes poor. If rapid- 
ity gap survival probabilities expected for point-like resolved photons are applied instead [24J, 
the overall normalisation is acceptable and the dependence of the data is better described. 
However, the description of the x J 7 ets dependence remains problematic. 

The analysis of the double-differential cross section d 2 cr/ dE^ n dx^ ts indicates that the E^ 1 
dependence of the data/theory ratio originates from the resolved photon-enriched region of x^ ts . 
However, the data are also consistent with no dependence on E^ n for either of the s J 7 cts re- 
gions studied. The ratio of the data to the NLO theory for the double-differential cross section 
d 2 cr/dz^, s dx J 7 ets is constant within errors throughout the region studied, indicating that the gap 
survival probability is insensitive to the presence or nature of remnants of either the photon or 
the diffractive exchange. 

Measurements of the ratio of diffractive to inclusive single-differential cross sections are 
presented as a function of several variables. The influence of multiple interaction effects in the 
inclusive data is large in the kinematic range studied here. The large uncertainties in modelling 
these multiple interactions preclude strong conclusions about rapidity gap survival on the basis 
of these data, although a reasonable description of the ratios can be obtained with suitably tuned 
Monte Carlo models. 
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Figure 2: Diffractive dijet photoproduction cross sections differential in (a) x 1 ' , (b) E 3 T , (c) 
logxjp, (d) z^ s , (e) <V cts ), (f) | A^ cts |, (g) W, (h) M 12 and (i) M x . The data points are shown 
with inner error bars corresponding to statistical uncertainties and outer error bars representing 
statistical and uncorrelated systematic uncertainties added in quadrature. The correlated sys- 
tematic errors are indicated by the open bands between the two solid black lines. The white 
lines show NLO QCD calculations obtained using the FR framework il26l[32ll33l and the HI 
2006 Fit B DPDFs, corrected for hadronisation effects. The dark bands around the theoretical 
predictions indicate the result of propagating the uncertainties on the Fit B DPDFs to the NLO 
calculation. The light bands show this DPDF uncertainty added in quadrature with the effect on 
the calculation of varying jiR and jip by factors of 0.5 and 2.0. In all figures, the predictions of 
the RAPGAP MC model are also shown. 
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Figure 3: Ratios of the single-differential cross sections to the corresponding NLO QCD predic- 
tions based on the FR framework and the HI 2006 Fit B DPDF set, corrected for hadronisation 
effects. See the caption of figure [2]for further details. 
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Figure 5: Single-differential diffractive dijet photoproduction cross sections as in figure [2l The 
FR theoretical prediction for resolved photons is modified by applying the scale factors from the 
KKMR model for point-like interactions ('KKMR suppressed') ll24ll or for hadron-like interac- 
tions ('resolved x0.34') [|23ll . The direct photon contribution is left unchanged in both cases. 
See the caption of figure [2] for further details. 
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Figure 6: Ratios of the single-differential cross sections to the corresponding NLO QCD pre- 
dictions based on the FR framework and the HI 2006 Fit B DPDF set, corrected for hadro- 
nisation effects. The FR theoretical prediction for resolved photons is modified by applying 
the scale factors from the KKMR model for point-like interactions ('KKMR suppressed') ll24l 
or for hadron-like interactions ('resolved x0.34') ll23l . The direct photon contribution is left 
unchanged in both cases. See the caption of figure [2] for further details. 
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Figure 8: (a)-(b) Double-differential cross section d 2 a/d^ s dx J 7 ets . (c)-(d) Ratio of the mea- 
sured double-differential cross section to the NLO prediction using the HI 2006 Fit B and HI 
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Table 2: Bin averaged hadron level differential cross sections for diffractive dijet photopro- 
duction. For each data point, the absolute statistical, uncorrelated and correlated systematic 
uncertainties and the hadronisation correction factors applied to the NLO calculations are given 

as S sm , (5 uncorr and S cori and 1 + (5 hadr , respectively. 
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Table 3: Ratios of differential cross sections for diffractive dijet photoproduction to NLO QCD 
calculations obtained in the FR framework with the HI 2006 Fit B DPDFs, corrected for hadro- 
nisation. For each ratio point, the absolute statistical, uncorrelated and correlated systematic 
uncertainties from the data and the absolute DPDF and scale uncertainties from the theory are 
given as 5 sm , 5 m , 5 corr , <W and 5 xa , respectively. 
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Table 4: Bin averaged hadron level double-differential cross sections for diffractive dijet photo- 
production. For each data point, the absolute statistical, uncorrelated and correlated systematic 
uncertainties and the hadronisation correction factors applied to the NLO calculations are given 
as 5 sm , 5 uncorr and 5 coir and 1 + <5 hadr , respectively. 
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Table 5: Ratios of double-differential cross sections for diffractive dijet photoproduction to 
NLO QCD calcualtions obtained in the FR framework with the HI 2006 Fit B DPDFs, corrected 
for hadronisation. For each ratio point, the absolute statistical, uncorrelated and correlated 
systematic uncertainties from the data and the absolute DPDF and scale uncertainties from the 
theory are given as S sUt , S mcon , <5 corr , 5 DPDF and 5 sca „ respectively. 
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Table 6: Ratios of the diffractive to the inclusive single-differential hadron level cross sections. 
The corresponding statistical, uncorrelated and correlated systematic uncertainties propagated 
to the ratio are given by 5 M , S mam and S mn , respectively. 
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